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Using the maximum entropy method, spectral functions of the pseudo-scalar and vector mesons are extracted 
from lattice Monte Carlo data of the imaginary time Green's functions. The resonance and continuum structures 
as well as the ground state peaks are successfully obtained. Error analysis of the resultant spectral functions is 
also given on the basis of the Bayes probability theory. 



1. Introduction 

The spectral functions (SPFs) of hadrons play 
a special role in physical observables in QCD (See 
the examples in However, the lattice QCD 

simulations so far have difficulties in accessing 
the dynamical quantities in the Minkowski space, 
because measurements on the lattice can only 
be carried out for discrete points in imaginary 
time. The analytic continuation from the imagi- 
nary time to the real time using the noisy lattice 
data is highly non-trivial and is even classified as 
an ill-posed problem. 

Recently we made a first serious attempt to 
extract SPFs of hadrons from lattice QCD data 
without making a priori assumptions on the spec- 
tral shape ||. We use the maximum entropy 
method (MEM), which has been successfully ap- 
plied for similar problems in quantum Monte 
Carlo simulations in condensed matter physics, 
image reconstruction in crystallography and as- 
trophysics, and so forth [|Jj|]. fn this report, we 
present the results for the pseudo-scalar (PS) and 
vector (V) channels at T = using the continuum 
kernel and the lattice kernel of the integral trans- 
form. The latter analysis has not been reported 
in§. 

2. Basic idea of MEM 

The Euclidean correlation function D(t) of an 
operator C(r, x) and its spectral decomposition 



at zero three-momentum read 



D(r) 



(O t (r,x)O(0,0))d 3 a: 

OO 

K(t, Lo)A(uj)du, 



(1) 



*Talk given by Y. Nakahara and M. Asakawa at 
LATTICE99. 



where r > 0, uj is a real frequency, and A(lo) 
is SPF (or sometimes called the image), which 
is positive semi-definite. The kernel K{t, lo) is 
proportional to the Fourier transform of a free 
boson propagator with mass u: At T = in the 
continuum limit, K = K cont (T,Lu) — exp(— tlj). 

Monte Carlo simulation provides D{ti) on the 
discrete set of temporal points < n/a < N T . 
From this data with statistical noise, we need 
to reconstruct the spectral function A(lu) with 
continuous variable to. This is a typical ill- 
posed problem, where the number of data is much 
smaller than the number of degrees of freedom to 
be reconstructed. This makes the standard like- 
lihood analysis and its variants inapplicable Q 
unless strong assumptions on the spectral shape 
are made. MEM is a method to circumvent this 
difficulty through Bayesian statistical inference of 
the most probable image together with its relia- 
bility |. 

MEM is based on the Bayes' theorem in prob- 
ability theory: P[X\Y] = P[Y\X]P[X]/P[Y], 
where P[X|Y] is the conditional probability of 
X given Y. The most probable image A{u>) for 
given lattice data D is obtained by maximizing 
the conditional probability P[A\DH], where H 
summarizes all the definitions and prior knowl- 
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edge such as A(uo) > 0. By the Bayes' theorem, 



P[A\DH] cx P[D\AH]P[A\H], 



(2) 



where P[D\AH] (P[A\H}) is called the likelihood 
function (the prior probability). 

For the likelihood function, the standard % 2 is 
adopted, namely P[D|yli/] = Z^ 1 exp(— L) with 

L = l^(D(n)-D A (r l )) (3) 



2 ^ x 



x Cr 1 1 (D(r j )-D A (r j )). 



Zl is a normalization factor given by Zl = 
(2 7 r) Ar /VdetC with N = r max /a - r mm /a + 1. 
D{n) is the lattice data averaged over gauge con- 
figurations and D (ri) is the correlation function 
defined by the right hand side of ([!]). C is an 
N x N covariance matrix of the data with N be- 
ing the number of temporal points to be used in 
the MEM analysis. The lattice data have gener- 
ally strong correlations among different t's, and 
it is essential to take into account the off-diagonal 
components of C. 

Axiomatic construction as well as intuitive 
"monkey argument" show that, for positive 
distributions such as SPF, the prior probabil- 
ity can be written with parameters a and m 
as P[A\Ham] — Zg 1 exp(aS). Here S is the 
Shannon-Jaynes entropy, 



A(u) - m(u)) - A(ui) k>£ 



(A(l 



(4) 
doj. 



\m(uj) 

Zs is a normalization factor: Zs = / e aS T>A. a 
is a real and positive parameter and m(u) is a 
real function called the default model. 

In the state-of-art MEM jj], the output image 
A ou t is given by a weighted average over A and 
a: 



Aout M 



A{lo) P[A\DHam}P[a\DHm] VA da 



A a {uj) P[a\DHm] da. 



(5) 



Here A a (uj) is obtained by maximizing the "free- 
energy" 

Q = aS-L, (6) 



for a given a. Here we assumed that P[A\DHam] 
is sharply peaked around A a (u>). a dictates the 
relative weight of the entropy S (which tends to 
fit A to the default model m) and the likelihood 
function L (which tends to fit A to the lattice 
data). Note, however, that a appears only in the 
intermediate step and is integrated out in the final 
result. Our lattice data show that the weight fac- 
tor P[a\DHm], which is calculable using Q 
is highly peaked around its maximum a = a. 
We have also studied the stability of the A out (ui) 
against a reasonable variation of m(uj). 

The non-trivial part of the MEM analysis is 
to find a global maximum of Q in the functional 
space of A(uj), which has typically 750 degrees of 
freedom in our case. We have utilized the singular 
value decomposition (SVD) of the kernel to define 
the search direction in this functional space. The 
method works successfully to find the global max- 
imum within reasonable iteration steps. 

3. MEM with mock data 

To check our MEM code and to see the de- 
pendence of the MEM image on the quality of 
the data, we made the following test using mock 
data, (i) We start with an input image Ai n (u>) = 
oj 2 pi n (oj) in the p- meson channel which simu- 
lates the experimental e + e" cross section. Then 
we calculate D in (T) from A in (uj) using eq.(Q). 
(ii) By taking D in {Ti) at N discrete points and 
adding a Gaussian noise, we create a mock data 
D m ock{Ti)- The variance of the noise cr(ri) is 
given by <j(ri) — b x Di n (ri) x r^/a with a param- 
eter b, which controls the noise level ^ . (iii) We 
construct the output image A out {us) = u) 2 p ou t(uj) 
using MEM with D rnock (T min < n < T max ) and 
compare the result with Ai n (ui). In this test, we 
have assumed that C is diagonal for simplicity. 

In Fig.l, we show pi n (u>), and p ou t{u) for two 
sets of parameters, (I) and (II). As for m, we 
choose a form m(uj) — itiqw 2 with iuq — 0.027, 
which is motivated by the asymptotic behavior 
of A in perturbative QCD, A(u> > IGeV) = 
(1 / Air 2 ) (1 + a s I 'it)uj 2 . The final result is, however, 
insensitive to the variation of uiq even by factor 5 
or 1/5. The calculation of A out {uj) has been done 
by discretizing the w-space with an equal sepa- 
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ration of 10 MeV between adjacent points. This 
number is chosen for the reason we shall discuss 
below. The comparison of the dashed line (set (I)) 
and the dash-dotted line (set (II)) shows that in- 
creasing Tmax and reducing the noise level b lead 
to better SPFs closer to the input SPF. 

We have also checked that MEM can nicely re- 
produce other forms of the mock SPFs. In par- 
ticular, MEM works very well to reproduce not 
only the broad structure but also the sharp peaks 
close to the delta-function as far as the noise level 
is sufficiently small. 
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Figure 1. The solid line is pi n (u>). The 
dashed line and dash-dotted line are p ou t(to) ob- 
tained with parameter set (I) a — 0.0847 fm, 
1 < r/a < 12, b = 0.001 and set (II) a = 0.0847 
fm, 1 < r/a < 36, b = 0.0001, respectively. 



4. MEM with lattice data 

To apply MEM to actual lattice data, quenched 
lattice QCD simulations have been done with the 
plaquette gluon action and the Wilson quark ac- 
tion by the open MILC code with minor modifica- 
tions §. The lattice size is 20 3 x 24 with (3 = 6.0, 
which corresponds to a = 0.0847 fm (a -1 = 2.33 
GeV), k c = 0.1571 jl0|], and the spatial size of the 
lattice L s a = 1.69 fm. Gauge configurations are 
generated by the heat-bath and over-relaxation 
algorithms with a ratio 1 : 4. Each configuration 



is separated by 1000 sweeps. Hopping parameters 
are chosen to be n = 0.153, 0.1545, and 0.1557 
with N con f = 161 for each n. For the quark prop- 
agator, the Dirichlet (periodic) boundary condi- 
tion is employed for the temporal (spatial) direc- 
tion. We have also done the simulation with pe- 
riodic boundary condition in the temporal direc- 
tion and obtained qualitatively the same results. 
To calculate the two-point correlation functions, 
we adopt a point-source at x — and a point-sink 
averaged over the spatial lattice-points. 

We use data at 1 < n/a < 12(24) for the 
Dirichlet (periodic) boundary condition in the 
temporal direction. To avoid the known patho- 
logical behavior of the eigenvalues of C Q], we 
take N conf > N. 

We define SPFs for the PS and V channels as 

A(u)=uj 2 p PS<v (u), (7) 

so that p ps v (u> — > large) approaches a finite 
constant as predicted by perturbative QCD. For 
the MEM analysis, we need to discretize the w- 
integration in (|l|). Since Aw (the mesh size) 
« 1/ r max should be satisfied to suppress the dis- 
cretization error, we take Aw = 10 MeV. uj max 
(the upper limit for the w integration) should be 
comparable to the maximum available momen- 
tum on the lattice: u> max ~ 7r/a ~ 7.3 GeV. 
We have checked that larger values of ui max do 
not change the result of A(ui) substantially, while 
smaller values of uj max distort the high energy end 
of the spectrum. The dimension of the image to 
be reconstructed is iV w = u> ma x/ Aw ~ 750, which 
is in fact much larger than the maximum number 
of Monte Carlo data N = 25. 

In Fig. 2 (a) and (b), we show the reconstructed 
images for each n in the case of the Dirichlet 
boundary condition. Here we use the continuum 
kernel K cont = exp(— tlu) in the Laplace trans- 
form. In these figures, we have used in = mow 2 
with m = 2.0(0.86) for PS (V) channel moti- 
vated by the perturbative estimate of mo (see 
eq.(^|) and the text below). We have checked 
that the result is not sensitive, within the sta- 
tistical significance of the image, to the variation 
of mo by factor 5 or 1/5. The obtained images 
have a common structure: the low-energy peaks 
corresponding to 7r and p, and the broad struc- 
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ture in the high-energy region. From the posi- 
tion of the pion peaks in Fig. 2(a), we extract 
k c = 0.1570(3), which is consistent with 0.1571 
fl(| determined from the asymptotic behavior of 
D(t). The mass of the p- meson in the chiral 
limit extracted from the peaks in Fig. 2(b) reads 
m p a — 0.348(15). This is also consistent with 
m p a = 0.331(22) jl0| determined by the asymp- 
totic behavior. Although our maximum value of 
the fitting range T max /d — 12 marginally covers 
the asymptotic limit in r, we can extract reason- 
able masses for tt and p. The width of n and p 
in Fig. 2 is an artifact due to the statistical errors 
of the lattice data. In fact, in the quenched ap- 
proximation, there is no room for the p-meson to 
decay into two pions. 

As for the second peaks in the PS and V 
channels, the error analysis discussed in Fig. 4 
shows that their spectral "shape" does not have 
much statistical significance, although the exis- 
tence of the non-vanishing spectral strength is 
significant. Under this reservation, we fit the po- 
sition of the second peaks and made linear ex- 
trapolation to the chiral limit with the results, 
m 2nd /m p = 1.88(8)(2.44(11)) for the PS (V) 
channel. These numbers should be compared 
with the experimental values: 1^^(1300) /m-p = 
1.68, and m p( i 4 50) / 'm p = 1.90 or m p ( 1700 )/m p = 
2.20. 

One should remark here that, in the standard 
two- mass fit of D(t), the mass of the second res- 
onance is highly sensitive to the lower limit of the 
fitting range, e.g., m 2nd /m p = 2.21(27) (1.58(26)) 
for r min /a — 8(9) in the V channel with f3 = 6.0 
|l(J . This is because the contamination from the 
short distance contributions from r < r m i„ is 
not under control in such an approach. On the 
other hand, MEM does not suffer from this diffi- 
culty and can utilize the full information down to 
Tmin / a = 1- Therefore, MEM opens a possibil- 
ity of systematic study of higher resonances with 
lattice QCD data. 

As for the third bumps in Fig. 2, the spectral 
"shape" is statistically not significant as is dis- 
cussed in Fig. 4, and they should rather be consid- 
ered a part of the perturbative continuum instead 
of a single resonance. Fig. 2 also shows that SPF 
decreases substantially above 6 GeV; MEM au- 
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Figure 2. Reconstructed image p ou t{w) for the 
PS (a) and V (b) channels. The solid, dashed, 
and dotted lines are for k — 0.1557, 0.1545, and 
0.153, respectively. For the PS (V) channel, mo 
is taken to be 2.0 (0.86). cu m ax is 7.5 GeV in this 
figure and Fig. 3. 



tomatically detects the existence of the momen- 
tum cutoff on the lattice ~ ir/a. It is expected 
that MEM with the data on finer lattices leads 
to larger ultraviolet cut-offs in the spectra. The 
height of the asymptotic form of the spectrum at 
high energy is estimated as 



p v (w~6GeV) 
4^2 ( 



(8) 



1 



n 



0.86. 



The first two factors are the qq continuum ex- 
pected from perturbative QCD. The third fac- 
tor contains the non-perturbative renormaliza- 
tion constant for the lattice composite operator. 
We adopt Z v = 0.57 determined from the two- 
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point functions at (3 — 6.0 |ll|] together with 
a s = 0.21 and k = 0.1557. Our estimate in 
eq.(^) is consistent with the high energy part of 
the spectrum in Fig. 2(b) after averaging over u>. 
We made a similar estimate for the PS channel 
using Z ps — 0.49 (l^] and obtained p ps (uj — 
6GeV) ~ 2.0. This is also consistent with Fig. 
2(a). We note here that an independent analy- 
sis of the imaginary time correlation functions [^) 
also shows that the lattice data at short distance 
is dominated by the perturbative continuum. 

In Fig. 3(a) and (b), the results using the lat- 
tice kernel Ki at are shown. Ki at is obtained from 
the free boson propagator on the lattice. It re- 
duces to K cont when a — > 0. The other parame- 
ters and boundary conditions are the same with 
Fig.2(a,b). The difference of Fig. 2 and Fig. 3 can 
be interpreted as a systematic error due to the 
hnitcncss of the lattice spacing a. 
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5. Error analsis 

The statistical significance of the reconstructed 
image can be studied by the following procedure 
Assuming that P[A\DHam] has a Gaussian 
distribution around the most probable image A, 
we estimate the error by the covariance of the 
image, — ((SaSaQ)~ 1 )a=A> where 5a is a func- 
tional derivative and (•) is an average over a given 
energy interval. The final error for A out is ob- 
tained by averaging the covariance over a with 
a weight factor P[a\DHm]. Shown in Fig. 4 is 
the MEM image in the V channel for n = 0.1557 
with errors obtained in the above procedure. The 
height of each horizontal bar is (/cw(w)) in each 
u> interval. The vertical bar indicates the error of 
(pout(ijj))- The small error for the lowest peak in 
Fig. 4 supports our identification of the peak with 
p. Although the existence of the non-vanishing 
spectral strength of the 2nd peak and 3rd bump 
is statistically significant, their spectral "shape" 
is either marginal or insignificant. Lattice data 
with better quality are called for to obtain better 
SPFs. 
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Figure 4. p ut{u) in the V channel for k 
0.1557 with error attached. 



Figure 3. Same with Fig. 2 except for the use of 
the lattice kernel K\ at . 
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6. Summary 

We have made a first serious attempt to recon- 
struct SPFs of hadrons from lattice QCD data. 
We have used MEM, which allows us to study 
SPFs without making a priori assumption on the 
spectral shape. The method works well for the 
mock data and actual lattice data. MEM pro- 
duces resonance and continuum-like structures in 
addition to the ground state peaks. The statisti- 
cal significance of the image can be also analyzed. 
Better data with finer and larger lattice will pro- 
duce better images with smaller errors, and our 
study is a first attempt towards this goal. 

There are many problems which can be ex- 
plored by MEM combined with lattice QCD data. 
Some of the applications in the baryon excited 
states, hadrons at finite temperature, and heavy 
quark systems will be reported in future publica- 
tions |3j. 
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